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Abstract. The instability in the cosmic-ray precursor of a supernova shock moving in interstellar medium is 
studied. The level of magnetohydrodynamic turbulence in this region determines the maximum energy of parti- 
cles accelerated by the diffusive shock acceleration mechanism. The high efficiency of cosmic ray acceleration is 
accepted, and the consideration is not limited by the case of weak turbulence. It is assumed that Kolmogorov 
type nonlinear wave interactions together with the ion-neutral collisions restrict the amplitude of random mag- 
netic field. As a result, the maximum energy of accelerated particles strongly depends on the age of supernova 
remnant. The maximum energy can be as high as ~ 10^^ Z eV in young supernova remnants and falls down to 
about ~ 10^'^ Z eV at the end of Sedov stage. Thus the standard estimate of maximum particle energy based on 
the Bohm limit calculated for the interstellar magnetic field strength is not justified in this case. This finding may 
explain why the supernova remnants with the age more than a few thousand years are not prominent sources of 
very high energy gamma-rays. 
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1. Introduction 

, Particle acceleration by supernova blast waves is com- 
' monly accepted as the main mechanism of acceleration 
. of galactic cosmic rays, see Drury et al. (pOOl), Malkov 
' & Drury H2U(J1|I for review. The high density of energetic 
, particles in the shock vicinity where diffusive shock ac- 
celeration occurs leads to the modification of gas flow 
through the action of energetic particle pressure and ini- 
tiates plasma instabilities produced by the current of en- 
ergetic particles. The plasma instabilities tend to increase 
the level of magnetohydrodynamic turbulence which scat- 
ter the particles. In the case of strong turbulence 6B ~ Bq 
((5B is the amplitude of random magnetic field, Bq is the 
regular magnetic field upstream of the shock), it may 
decrease the cosmic ray diffusion coefficient D to the 
Bohm value Db = vrg/i, where v is the particle velocity, 
Tg = pc/^ZeBo) is the Larmor radius of a particle with 
momentum p and charge Ze. The Bohm value of diffusion 
coefficient is commonly used in simulations of cosmic ray 
acceleration in supernova remnants, e.g. Berezhko et al. 
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The dependence of diffusion on energy determines 
the maximum energy which particles can gain in the 
process of acceleration. The necessary condition of ef- 
ficient acceleration at the shock at the free expansion 
stage is D{E) < O.lushRsh, where Rsh is the radius and 
Ugh — dR/dt is the velocity of spherical shock. Notice that 
the typical supernova burst with kinetic energy of ejecta 
£ = SqiIO^^ erg in the interstellar gas with number den- 
sity riQ cm^'^ gives the maximum value of the product 
O.lUshRsh ~ 10^''(f5i/no)° '* cm^ s~^ at the end of the 
free expansion stage of supernova remnant evolution when 
this product reaches its maximum value. At the same time, 
the typical value of cosmic ray diffusion coefRcient in the 
Galaxy is Dism = 5.9x lO^^/Ji?^^ cm^ s~\ where /3 = v/c 
is the dimensionless particle velocity and Rm — pc/Z is 
the particle magnetic rigidity in units GV, see Jones et al. 
(2001), and thus the necessary condition of acceleration 
can not be fulfilled for relativistic particles unless their 
diffusion coefRcient is anomalously small in the vicinity 
of the shock. The Bohm value which is a lower bound of 
diffusion coefficient along the magnetic field at SB < Bq 
is equal to Db = 6 x 10'^^ PR^ cm'^/s at = 5 x 10"^ 
G. The entire pattern of cosmic ray acceleration on su- 
pernova shocks critically depends on the assumption that 
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the energetic particles themselves produce the turbulence 
needed to provide the anomalously slow diffusion at the 
cite of acceleration. This was perceived at the very begin- 
ning of the studies of diffusive shock acceleration mecha- 
nism rBell. [W751 Lagage & Cesarsky,[TO83). 

Analyzing the early stage of supernova remnant evo- 
lution when the shock velocity is high, Ush ~ 10* km s""'^, 
Lucek & Bell f2000,), Bell & Lucek l^MH^ found that the 
cosmic ray streaming instability can be so strong that the 
amplified field SB > 10~* G far exceeds the magnetic field 
in the undisturbed medium ahead of the moving shock 
Bq ~ (2-f-5) X 10~^ G (see also earlier papers of McKenzie 
& Volk itim . Vol k et al . (HM . Falle & Giddings (tlM . 
Bennett & Ellison H1995|l 'l. The maximum particle energy 
increases accordingly. The cosmic ray acceleration contin- 
ues at later stage of supernova evolution when the mass 
of swept-up interstellar gas exceeds the mass of super- 
nova ejecta and the shock velocity goes down according 
to Sedov solution. The major portion of galactic cosmic 
rays is accelerated by the expanding blast wave during 
this stage. The cosmic-ray streaming instability is less ef- 
ficient as the shock velocity falls down with time and the 
nonlinear wave interactions can significantly reduce the 
level of turbulence at late Sedov stage (Volk et al. 1988 
Fedorenko ll99()|l . This leads to fast diffusion and decreases 
the maximum energy which can be achieved in a process 
of diffusive shock acceleration. This effect is aggravated 
by the wave damping which is due to the ion-neutral col- 
lisions if neutral atoms are present in the interstellar gas 
upstream of the shock fBell. IW^ Drury et al. [TMH| . 

In the present paper, we consider the acceleration of 
cosmic rays and their streaming instability in a wide range 
of shock velocities. The level of magnetic field fluctuations 
is allowed to be arbitrary large {SB Bq), and the rate of 
nonlinear wave interactions is assumed to correspond the 
Kolmogorov type nonlinearity. The dissipation of waves 
by the ion-neutral collisions is also taken into account in 
the case of supernova burst in not fully ionized interstellar 
gas. The task is to find the maximum energy of acceler- 
ated particles as a function of age of supernova remnant. 
The problem is enormously complicated by the fact that 
the processes of wave-particle and wave- wave interactions 
have no universally accepted theoretical description in the 
case of strong turbulence. So, our consideration is made 
at the level of approximate estimates and we present not 
more than a schematic description of the overall picture. 



2. Maximum energy of accelerated particles 

We consider the adiabatic stage of supernova remnant ex- 
pansion when the major fraction of cosmic rays injected 
in the interstellar medium is accelerated. The evolution 
of a shock produced by supernova outburst in interstellar 
medium is described by Sedov solution, see e.g. Lozinskaya 
(UnSSl, Ostriker & McKee ifTU^ . In the case of very effi- 
cient acceleration when the cosmic ray pressure inside the 
supernova envelope governs its expansion, the equations 



for Rshit) and Ush{t) are: 



Rsh = 4.26 



no 



VlOV 



2/5 



pc, 



Ush = 1.67 X 10'' 



£51 
no 



1/5 



103yr 



-3/5 



kms-^ (1) 



The selfsimilar solution (1) is valid under the con- 
dition that the mass of swept-up gas significantly ex- 
ceeds the mass of ejecta Alej (^ ^Mq) that is fulfilled 



at Rsh > Ro 



(3Mej/47rm^no)^/^ 



1.9n, 



-1/3 



pc 



t > tQ — R/uq = 190riQ yr, where uq is the num- 
ber density of hydrogen atoms in the interstellar gas, 
nia = 1.4m is the mean interstellar atom mass per hy- 
drogen nucleus, m is the proton mass, and uq ^ 10^ km 
s^^ is the ejecta initial velocity. The preceding evolution 
of a supernova remnant can be approximated as a free 
expansion with velocity uq. The adiabatic approximation 
used in the derivation of equations (1) is not valid when 
the shock velocity sinks below about Ugh = 280£^f*^^nP'^^^ 
km and the radiation energy losses become essential. 
The evolution of a supernova shell with intense radiative 
coofing goes as Rsh oc t^/'', Ush oc t~^/'^ . We set £"51 = 1 
in the subsequent numerical estimates. 

In the following, we distinguish tiq, the mean gas num- 
ber density which determines the shock wave expansion 
law, and n, the intercloud gas number density. It corre- 
sponds to cloudy structure of interstellar medium with 
small dense clouds randomly distributed in the diffuse 
background gas. Two phases of interstellar gas are con- 
sidered below: a hot ionized component with the temper- 
ature of background gas T = 7 x lO^K, the number den- 



sities no = 10 cm 



3 X 10 cm and the mag- 



netic field strength Bq = 2 /iG; and a warm partly ionized 
component with the following parameters: T = 8 x lO^K, 
no = 0.4 cm~^, n = 0.1 cm^^, n^ — 0.03 cm~^, Bq — 5 
fiG (where n,j is the number density of ions). We assume 
that the diffusive shock acceleration of energetic particles 
should cease to operate after 10^ years in the hot gas when 
the shock Mach number falls below M ~ 3, and by the 
time 5 x 10* yr in the warm gas when Sedov stage ends. 

The acceleration of a fast particle diffusing near the 
shock front is a version of Fermi type acceleration. The 
fast particles are scattered on inhomogeneities of magnetic 
field frozen into background plasma and gain energy cross- 
ing the shock where plasma is compressing. In the test par- 
ticle approximation, the distribution of accelerated cosmic 
ray particles on momentum for high Mach number shocks 
has the canonical power- law form f{p) oc up to some 
maximum momentum Pmax- In the case of efficient accel- 
eration, the action of cosmic ray pressure on the shock 
structure causes nonlinear modification of the shock that 
changes the shape of particle spectrum making it fiatter at 
relativistic energies, see Berezhko & Ellison (1999). With 
allowance made for the nonlinear modification, we assume 
in the following that the spectrum of energetic particles at 
the shock is of the form f{p) cc p~*+'^ where < a < 0.5. 
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The normalization of function f{p) is such that the inte- 
gral N = An J dpp^f{p) gives the number density of rela- 
tivistic particles. We assume that the pressure of cosmic 
rays at the shock P = (47r/3) / dpp^vf{p) is some frac- 
tion < 1 of the upstream momentum flux entering the 
shock front, so that P = S,crpu1f^ {p = maUo). The typ- 
ical value of ^cr ~ 0.5 and the total compression ratio 
7 were found in the numerical simulations of strongly 
modified shocks by Berezhko et al. (1996) under the as- 
sumption of efficient dissipation of Alfvenic turbulence in 
the shock precursor. Thus we accept the following distri- 
bution function of cosmic rays at the shock: 



foip) = 

<<5(Pmax) = 



'i^crPult^H{pn^^y, ~p) 1 

47rC(/9(pmax) [mcYp-^-" ' 



(2) 



where H{p) is the step function. The ultra relativis- 
tic asymptotics of Lp{p) is •^(jp) ~ {p/mcY/a and 
the good approximation at p/mc 1 is '/'(Pmax) ~ 
a"HPmax/mc)"(l-((l+a)(p/mc)'')-i). The value a = 0.3 
is used in the numerical estimates below. 

The spatial distribution of accelerated particles in the 
plane shock approximation can be roughly presented as 



h{p,x) = /o(p)exp Ush 



D{p,xi) 



at X < (upstream of the shock); and 
f2{p,x) = fa{p) 



(3) 



(4) 



at X > (downstream of the shock). 

Here we use the reference frame where the shock is at 
rest at X = 0, the gas velocity is parallel to the x axis, and 
Ux > 0. The steady state one-dimcnsional solution (3) is 

1/2 

valid for |x| < mm{Rsh, Rsh {D/uR) ' }. Equation (3) 
shows the presence of cosmic ray precursor of characteris- 
tic size D/ush ahead of the shock. 

The cosmic ray diffusion coefficient is determined by 
the following approximate equation which can be used for 
an arbitrary strong turbulence, 5B ^ B^: 

,1/2 



D 



3^2 (> kr 



~^^g-! ^res 



l + Ai 



(5) 



We introduce here the dimensionless amplitude of the 
random field A^^ot = SB^/B^ = (47r/B2) / dkW{k) , 
where W{k) is the energy density of magnetohydrody- 
namic turbulence. The Larmor radius is defined through 
the regular field rg ~ pc/ZeBo- The scattering of par- 
ticles of Larmor radius rg is mainly produced by ran- 
dom inhomogcneities with the resonant wave number fc^es 
and the amplitude of resonant waves is characterized by 
A'^{> kres) — 4:TTkresW{kres)/BQ. Wc Set up the formula 
(5) with the rounded numerical coefficient as the gener- 
alization of equation for cosmic ray diffusion coefficient 



in weak random field (Berezinskii et al.. I199(j|l . In doing 
the generalization, we distinguish for each energetic par- 
ticle with momentum p a large scale magnetic field, which 
includes the regular magnetic field plus the averaged ran- 
dom magnetic field with k < kres{p), and a small scale 
magnetic field, which includes the random magnetic field 
with k > kres{p)- The particle moves adiabatically in a 
large scale field and experiences scattering by small scale 
field. Note that the choice of numerical coefficient in (5) 
affects the numerical coefficient in the equation (7) below. 

It is usually assumed that the random field which is 
necessary for cosmic ray diffusion inside the supernova 
shell can be produced by the advection of turbulence 
through the shock front, by the intrinsic instabilities at 
the shock transition region, and by the irregular gas flow 
downstream of the shock. The turbulence in the upstream 
region should be due to the cosmic ray streaming insta- 
bility. 

In the case under study, the following steady-state 
equation determines the energy density W{k, x) of the tur- 
bulence amplified by the streaming instability in cosmic 
ray precursor upstream of the shock: 



rni)w. 



(6) 



X < 0. 

Here the left hand side of the equation describes the 
advection of turbulence by gas flow with highly supersonic 
velocity u = Usu- The terms on the right hand side of the 
equation describe respectively the wave amplification by 
cosmic rays, the linear damping of waves in background 
plasma, and the nonlinear wave-wave interactions which 
limit the amplitude of turbulence. In the strict sense, r„; is 
some integral-differential operator acting on the function 
W{k,x). 

The diffusion current which is due to the nonuniform 
distribution of cosmic ray particles upstream of the shock 
causes the streaming instability and the amplification of 
magnetohydrodynamic waves with the growth rate 



c^k 



dpp 1 



Pres {k) 



D 



Pr.s(k) 



df 



dx 



(7) 



(the waves with kx < ^ are amphfied). Here Pres{k) = 

1 /2 

ZeSo (l + Afot) /ck, and the Alfven velocity is de- 
fined through the regular field Bq as Va — Bq / y/Anp = 
l.S{BQ/l^G)n^^^^ km s^^. The most abundant proton 
component of cosmic rays determines the value of Per, so 
we put Z = I in the subsequent calculations. The expres- 
sion (7) is the generalization of equation for weak random 
field (Berezinskii et al, 1990) to the case of strong random 
field. 
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Equations (2), (3), (5), (7) let us estimate the value of 
growth rate in the upstream region {x < 0) as 



CM^.rUl^k'-'^ (1 + Al,) 

kA'^{> k,xi) 



-(l-a)/2 



exp 3ush / dxi 



(8) 



at A; > fcres(Pmax)j and = Otherwise. Here rgo 



27 



4(5-a)(2-a) 



^, SO that Ccr(0.3) 



mc^/{eBo), Ccr{a) 
0.845. 

Considering that all waves upstream of the shock are 
generated by cosmic-ray streaming instability, the parti- 
cles with maximum momentum Pmax are in resonance with 
the waves which have the minimum wave number k^i^i = 
(1 + AltY^^ r-i(Pmax). Note that A^{> k^ir.) = ^L- 

The most important mechanism of linear damping of 
Alfvenic turbulence in the problem under consideration is 
the wave dissipation due to ion-neutral collisions, see e.g. 
Kulsrud & Cesarsky (1971). It occurs in not fully ionized 
gas and is described by the equations 



(9) 



r, - 



fc2(l 



^tot 



2i^i„(l -I- ni/nny 



AltVa 



where = nn (vtha) « 8.4 x IQ-^ (T/IO^K)^' (jih/I 
cm~^) s~^ for temperatures T ^ lO^K to lO^K is the 
frequency of ion-neutral collisions with the cross section 
a averaged over velocity distribution of thermal particles, 
riH is the number density of neutral hydrogen. This gives 
Vin — 7.7 X 10^^° s^^ in warm diffuse gas. 

In spite of a great progress in the investigation of mag- 
netohydrodynamic turbulence, it is not easy to specify Tni 
in equation (6). The incompressible MHD simulations of 
driven and decaying turbulence by Verma et al. H1996|l 
demonstrated the presence of Kolmogorov type energy 
cascade to large wave numbers that can be qualitatively 
described by the equation 



TnlW^ik) 



dk 



1/2 



W^{k) 



(10) 



where Ck ~ 3.6 is the so-called "Kolmogorov constant" 
and ± signs correspond to Alfven waves travelling in 
opposite directions, i.e. with k^ > and fc^, < re- 
spectively, the total wave energy density is W{k) — 
W^{k) + W^{k), and the turbulence is assumed to 
be weak, kW^{k) (i3o/47r) ^ 1. The nonlinear cas- 
cade of Alfvenic waves is anisotropic. The main part 
of energy density in this turbulence is concentrated at 
perpendicular to the local magnetic field wave vectors 
fc_L k, while the parallel wave numbers are small: 



1/2 



k±, see Goldreich & Sridhar 

l(TM5|l . Cho & Vishniak iPTOjl . 

The theoretical approach based on the kinetic theory of 
weakly turbulent coUisionless plasma gives somewhat dif- 
ferent results, see Livshits & Tsytovich fl970"), Lee & Volk 
|1973), Kulsrud (1978), Achterbcrg (1981), Fedorenko et 
al. (1990), Zirakashvili (2000|. In particular, it was shown 
that the nonlinear interactions of Alfven waves and ther- 
mal particles lead to the nonlinear wave damping which 
is described by the equation 



TrdW^ik) 



f kW^jk) 



W^{k) 



(11) 



in the case of high f3th = v'^/Va > 1 plasma (Zirakashvili, 
I2()()()|l and by the same within an order of magnitude 
expression but without factor \/ (3th in the case of low 
Pth < 1 (Livshits & Tsytovich, 1970 ZirakashviH, ESni). 
Equation (11) was obtained for the case of weak turbu- 
lence and under the assumption that wave distribution is 
not strongly anisotropic in fc-space. A rate similar to (11) 
was derived for the wave-wave interactions in magneto- 
hydrodynamics based on the theory of perturbations at 
small wave amplitude (Chin & Wentzel. 119721 Akhiezer et 
aL.IW^ SkiUing.HTO. 

The observations show that the turbulence in the solar 
wind plasma is of Kolmogorov type, that evidently favors 
equation (10) to equation (11), and contains Alfven waves 
moving both along the magnetic field and in perpendicular 
directions (Saur & Bieber,'199Q|. The data on interstellar 
turbulence are also consistent with the assumption that 
a single Kolmogorov spectrum extends from scales 10* to 
10^" cm (Armstrong et al. . 1199511 . 

It is clear from expressions (10) and (11) that Alfven 
waves with the opposite signs of kx are needed for 
nonlinear interactions. (The only exception is the one- 
dimensional case when the interactions are allowed for 
the waves propagating in one direction along the external 
magnetic field, see Lee & Volk ((TF73), Kulsrud lfT?T7^ . 
Achterberg l|1981|l and Zirakashvili l|2000,) for correspond- 
ing equations and discussion.) The cosmic-ray streaming 
upstream of the shock amplifies waves W~ with fc^, < 
that propagate in the direction of cosmic ray diffusion flux, 
and damps the waves with k^ > 0. So, one needs an 
additional mechanism to maintain the waves in both direc- 
tions and to assure their efficient nonlinear interactions. 
This is provided by the scattering of waves on the fluc- 
tuations of gas density. Considering the case of the most 
strong wave interactions, we assume in the following that 
the isotropization of waves is fast and the waves mov- 
ing in both directions along the external magnetic field 
are present. A similar approach is used in the studies of 
Alfven waves in the solar wind turbulence, see Goldstein 

mm. 



et al. fT995), Hu et al. 

The Kolmogorov type nonlinearity is assumed in the 
consideration that follows. The simplified expression 



r„, = {2CKy''/^VakA{> k) « 0.05Kfc^(> fc), 



(12) 
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at Ck = 3.6 is accepted in our calculations. The order of 
magnitude estimate (12) is based on equation (10) and is 
valid for an arbitrary strong turbulence. 

It is worth noting that the earlier consideration of 
cosmic ray streaming instability in supernova remnants 
by Lagage & Cesarsky (HHHSll, Volk et al. (1988) and 
Fedorenko l|199()(l was made in the limit of 5B < Bq and 
with the use of equation similar to (11) that is very dif- 
ferent from our present approach. 

Let us consider separately the cases of nonlinear and 
linear dissipation of waves. 

A. Cosmic-ray streaming instability with Kolmogorov 
type nonlinearity. 

The linear damping is ignored in this case. The wave 
density is determined by equation (6) at F; = and with 
Vni given by equation (12). The exponential distribution 
of cosmic rays (3) with a characteristic scale D/ush de- 
termines the spatial distribution of waves in the upstream 
region that gives the estimate of the first term of equa- 
tion (6) at the shock udW/dx « u^W/D in the case of 
efficient wave amplification. Now using equations (5), (8), 
and (12), one can find from (6) the following approximate 
equation that determines the amplitude of random field 
at the shock: 



cVaf{Pmax){krgo)'' (1 



Alt 



,(l-i)/2- 



(13) 



Equation (13) allows to find the following equation for 



A, 



'^11'^ /42 



2v{l 



+ {2CKY^I^VaAtot 



cVaV (Pmax)(p 11 



,1/2- 



(14) 



As was pointed out in the Introduction, the efficient 
particle acceleration at the shock necessitates sufhciently 
small diffusion coefficient. The maximum momentum of 
accelerated particles at the Sedov stage of SNR evolution 
is approximately determined by the condition -D(pmax) — 
KUshRsh, where x « 0.04, see e.g. Berezhko ct al. (1996). 
The particles with larger diffusion coefficient are not accel- 
erated and not confined in the remnant. The last condition 
and equation (5) give 



Pn 



3KAlfUshRsh 



(1 



Al,) 



1/2, 



(15) 



rgo 



As the most abundant species, protons mainly drive 
the instability. Other particles experience resonant scat- 
tering on produced turbulence and undergo diffusive shock 
acceleration. As a result, the ions with charge Z have the 
same shape of spectra as protons but written as function 
oi p/Z instead of p. This is of course valid for the maxi- 
mum momentum attained at the shock and Pmax/Z should 



appear in equation (15) instead of Pmax if ^ 1 (as be- 
fore, m is the proton mass, rgo is the Larmor radius of a 
proton with momentum mc in equation (15)). The very 
high energy electrons are accelerated with the same rate 
as protons at the same magnetic rigidity. However, the 
maximum energy of electrons is smaller then for protons 
because of heavy energy losses on synchrotron radiation 
and inverse Compton scattering, see analysis of these pro- 
cesses by Gaisser et al. H1998|l . Reynolds H1998|l . 

Equations (14), (15) together with equations (1) de- 
termine the maximum momentum of accelerated particles 
(protons) during Sedov stage of SNR evolution in an im- 
plicit form. The results of numerical solution of these equa- 
tion for Pmax as the function of shock velocity Ush are pre- 
sented in Figures 1 and 2 for the cases of supernova rem- 
nant expansion in hot and warm interstellar gas respec- 
tively. The values of ^cr = 0.5 , a = 0.3, and £51 = 1 were 
accepted in these calculations. For comparison, the solu- 
tion without wave damping, F; = 0, F„i = is presented. 
The Bohm limit for maximum momentum calculated at 
the interstellar magnetic field strength Bq is also shown 
in Figures 1, 2. The Bohm limit represents the standard 
estimate of maximum momentum of accelerated particles 
that is commonly accepted in the investigations of cosmic 
ray acceleration in supernova remnants. Our considera- 
tion shows that this " standard " estimate of Pmax might 
be completely unsuitable. It does not work in young su- 
pernova remnants at the age i ^ 3 x 10^ yr when cosmic 
ray stream instability generates strong random magnetic 
field SB > Bq that reduces particle diffusion (in agree- 
ment with Bell & Lucek (2001) results). It also does not 
work in old supernova remnants at i > 3 x 10'^ yr when 
strong Kolmogorov type nonlinear interactions of waves 
damp the turbulence generated by instability and thus in- 
creases particle diffusion upstream of the shock. 

It is useful to give simple analytic expressions for Pmax 
in the limits of high and low values of shock velocity Ush- 

It is easy to check that if the shock velocity is high 
enough so that 



Ush > TFTi'^sh > 



3Va 



9(2Ck)^/^ 



C 



(16) 



then the first term in the left-hand side of equation (14) 
dominates over the second term, and the wave amplitude 
is large. At » 1. Thus the advection of waves is more 
important than the nonlinear wave interactions and the 
turbulence is strong when (16) is fulfilled. The inequalities 
(16) are satisfied at Ug^ ^ 8.3 x 10^ cm s^"'^, Ug^ ^ 7.9 x 
10^ cm s""'^ in hot gas, and Ugh 3> 8.3 x 10^ cm s^^, Ug}^ 3> 
3.4 X 10^ cm s^^ in warm gas at a = 0.3, ^cr — 0.5. The 
equations (14), (15) now give 



Pn 



(17) 



mc rgoVaC 

2.2 X 105 iLr/0.5)£!i%'^'n^/'t^li'. 

Here we set x = 0.04, a — 0.3, Ck = 3.6, use equations 
(1) for Ush{t) and Rgh{t), and approximate (y3(pmax) « 
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(jPiaim/iTicY / a, V — c; ^Kyr IS the supernova remnant 
age in units 10'^ yr. Equation (17) estimates Pmax/wc « 



2 X lO^t^yl^ in hot gas, and Pmax/m 
warm gas at txyr ^ 3, £51 = l,^cr — 0.5. The cosmic 
ray diffusion coefficient in the upstream region depends 
on rigidity as £> oc r^^" at p < Pmax- 

In the opposite hmit of iow shock veiocity 



1 X lO^tKyr' in 



Ush < 



C 



1/5 



1/3 



(18) 



the second term in the ieft-hand side of equation (14) dom- 
inates over the first term and the wave amphtude is smail, 
At <^ 1. The non-iinear wave interactions are more impor- 
tant than advection in this case. The inequalities (18) are 
satisfied at <C 7.4 x 10'' cm s~^, Ugh 8.0 x 10'' cm s~^ 
in hot gas, and Ush <C 4.4 x lO"^ cm s~-^, u^h < 4.5 x 10^ 
cm s~^ in warm gas. The equations (14), (15) give 



0.2i^a^Cl{a)CU 



mc 



4.1 X 10^^ i^,r/0.5f£U'n-'^'n'B-U 



-19/5 



0,A'G''Kyr 



(19) 



7 X lOH^l^J^ in hot 



i/mc w 1 X lO^tj^yl'^"' in warm gas at txyr ^ 



that leads to the estimate p 
gas, and p, 

10. The cosmic ray diffusion coefficient in the upstream 
region depends on rigidity as D oc wr^"^" at p < Pmax in 
this case. 

As it is evident from Figures 1, 2 and equations (17), 
(19), the value of p,„ax is rapidly decreasing with the age 
of supernova remnant t. The maximum energy of cosmic 
rays accelerated by supernova shock at early Sedov stage 
is close to 3x lO^'' eV for protons and it falls down to about 
10^° eV at the end of Sedov stage for the set of parameters 
accepted in the present work. In particular, the particle 
energy is less than 10^^ eV for t > 5 x lO'^ yr and this may 
explain the absence of TeV gamma-ray signal from many 
supernova remnants (Buckley et al. I1998|l where gamma- 
rays could be produced through 7r° decays if sufficiently 
energetic protons and nuclei were present. The commonly 
accepted estimate of maximum particle energy based on 
the Bohm diffusion limit enormously overestimates Pmax 
for these objects if the level of self-generated turbulence is 
limited by the strong Kolmogorov type nonlinearity (12). 

B. Cosmic-ray streaming instability in warm gas with 
collisional wave dissipation. 

The cosmic ray acceleration in supernova remnants 
in the presence of wave dissipation by ion-neutral colli- 
sions was investigated by Drury et al. (1996). Different 
aspects of particle acceleration by supernova shocks mov- 
ing in partly ionized interstellar gas were also consid- 
ered by Boulares & Cox ifTUSSjl . Bykov et al. ^UUU^. It 
was assumed that the turbulence generated by cosmic ray 
streaming instability is weak. Here we do not restrict our 



consideration to the case of weak turbulence. We ignore 
the nonlinear wave dissipation in this section and consider 
the collisional dissipation in partly ionized warm gas. The 
wave density W{k) is determined now from equation (6) 
where the term r„; is omitted and the term F; is given by 
equation (9). So, similarly to the case A and associated 
equations (13), (14), we have the following equation for 
the amplitude of random magnetic field at the shock: 



3ul^AH> k)k 



(l-a)/2 



(20) 



cVa^{Pn,e.^){krgnY (1 + ^L) 
at fc > fc,, /C* = Vin/ ^LK) • 

Now, using equations (5), (15), (20) and relations 

A (> ^min) — A^Qt^ ^min 

following equation for pmax: 

2C,,(a)e cr ^sh 



^tot/^g,max 



we find the 



V mc / 



c) 



(21) 



vahd for p^ax < P*, p*/mc = (l 

The term i/i„ should be omitted in equation (21) if Pmax > 
. The dissipation on neutrals is not essential for the 
development of instability if the last inequality is fulfilled. 

The value of maximum particle momentum Pmax as 
a function of the supernova shock velocity Ugh deter- 
mined according to equation (21) is shown in Figure 2. 
The approximate expression based on equation (21) at 
lyin > 20ush/Rsh is 

Pma-K ^ 2(2Ccr(^)Ccr^s/i ^ 
mc cVargoVrn 

3 X 10^ (ec./0.5) sTn-'/'n-'^h^li' (22) 

at a = 0.3 that gives (pmax/'^ic) « 5.7 x 
10^ (it./j/lO^km s-i)^ « 1.64 X lO-^^Kyr^ at = 0.5, 
£51 = 1. The last approximate formula and Figure 2 ev- 
idence that the maximum energy of accelerated particles 
decreases from about 3 x 10^^ eV in early Sedov stage to 
3x 10^° eV at the end of Sedov stage due to the damping of 
cosmic-ray generated turbulence by ion-neutral collisions. 

The solution with a smaller value of Pmax at given 
shock velocity Ush should be chosen when both nonlinear 
and collisional wave damping are present. 

It is known that some envelopes of Type la supernovae 
expand in partly ionized interstellar gas (see Lozinskaya, 
1992), so that the equations obtained above are appropri- 
ate in this instance. On the other hand, an explosion of 
Type II supernova leads to the full ionization of surround- 
ing diffuse gas over distances of several of tens parsecs that 
holds about lO^n"-'^ years after the burst. Therefore, the 
intercloud neutral gas density is very low shortly after the 
burst and rises linearly in time because of recombination 
as 



riH = CrU t, 



(23) 
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Fig. 1. The maximum momentum of accelerated protons Pmax in units of mc as a function of shock velocity Ush at 
the Sedov stage of supernova remnant evolution in the hot interstellar gas is shown by solid line. The kinetic energy 
of the explosion is £ — 10^^ erg; the efficiency of cosmic ray acceleration ^cr — 0.5. The parameters of interstellar 
medium are the following: the gas temperature T = 7 x 10^ K, the average gas density no = 10~^ cm~^, the density 
of diffuse gas n = 3 x 10~^ cm"'^, and the interstellar magnetic field strength Bq = 2 /iG. The dash line presents the 
age of a supernova remnant t (ploted on the right ordinate) as a function of shock velocity. Also shown: dot line - the 
Bohm limit on maximum particle momentum in the field Bq; dash-dot line - the maximum particle momentum when 
the wave dissipation is not taken into account. 



at t <^ (Crn) , where Cr{T) is the recombination rate 
of hydrogen. At a later time t ~ {Crn)~ the value of 
hh tends to some constant value which is determined by 
the balance between ionization and recombination. The 
standard value of recombination rate Cr = 4.1 x 10~^^ 
cm-^s-i at T « lO^if (Seaton, TSS^ gives the time- 
dependent frequency of ion-neutral collisions Vin{t) — 
{vthfy) CrU^t « 1 X 10~-^°n^iKyr s~^. One can continue to 
use equations (20)- (21) with the substitution = Vinit) 
for the determination of maximum momentum of acceler- 
ated particles. Figure 2 illustrates the results. It is evident 
that the ionization of interstellar atoms by supernova ra- 
diation diminishes the coUisional wave dissipation. 



3. Supernova explosion in a stellar wind 

Immediately after the supernova burst, the shock propa- 
gates not in the interstellar medium but through the wind 
of a progenitor star. The duration of this phase depends 
on the wind power and it may cover even the considerable 
part of the adiabatic (Sedov) stage of supernova remnant 
evolution if the pre-supernova star was very massive, more 
than about 15 solar masses. The massive progenitors of 
Type-II and lb supernovae undergo a significant mass loss 



prior to explosion during their red and blue giant phases, 
and the Wolf-Rayet phase for the most massive stars. The 
stellar wind modifies the ambient density distribution and 
creates an expanding shell. The cosmic ray acceleration 
in the case of shock propagation through the wind was 
considered by Volk and Biermann Cl988| , Berezinsky and 
Ptuskin (1989), and Berezhko and Volk iPHnOjl under the 
assumption of Bohm diffusion for energetic particles. Bell 
and Lucek ()2()()l|l included the effect of strong cosmic-ray 
streaming instability into the consideration of high veloc- 
ity shocks. 

A distinctive property of the problem is the non- 
homogeneous distributions of gas density and magnetic 
field in the stellar wind fiow. In particular, the spheri- 
cally symmetric distribution of gas density in the stellar 
wind is n^,, = M/ (47rrnau™r2), where M = 10~^M_5 
Mq yr~^ is the mass loss, u^, = lO^Utu^g cm is the 
wind velocity, and the normalization and numerical val- 
ues here and below are given for the typical supernova 
Type II progenitor, the red giant star (see Berezinsky 
and Ptuskin, 1989). Similar to the interplanetary mag- 
netic field, the stellar wind magnetic field has the shape 
of a Parker spiral (Parker, '1958). The magnetic field has 
predominately azimuthal structure at distance from the 
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Iff Iff 10^ 

Uji^, km/s 

Fig. 2. The same as in Figure 1 but for the supernova burst in warm interstellar gas {T = 8 x 10"^ K, no — 0.4 cm"^, 
n — 0.1 cm~'^, Bq — 5 /iG). Three solid lines correspond to three cases of wave dissipation considered separately: 
nonlinear wave interactions; damping by ion- neutral collisions at given above number density of neutral atoms; damping 
by ion-neutral collisions when the diffuse neutral gas restors its' density after complete ionization by the radiation 
from supernova burst. 



star r u^/ {^Is'mO) « 5 x 10'^^/sin6' cm where its' value 
is Bq = B^,rlil sin 9 / (u^r). Here i?* = 10 G is the surface 
magnetic field strength at star radius — 2 x 10^^ cm , 
n — 2 X 10"^ s~^ is the angular velocity of star rotation, 
and 9 is the polar angle. Hence BQ(r)r « 8 x 10^'' sin 6* 
Gxcm at r ^ 5 X 10^^ cm. It is worth noting that the 

value of Alfven velocity 14 = BQ{r)r^J u^/M does not 
depend on distance r in this case. 

The shock propagation in a stellar wind obeys the fol- 
lowing equations (e.g. Lozinskava il992;) : 



R.h = 8.77 



Ush 



M-5 

5.72 X 10" 



1/3 



,2/3 



^(^)'''C^kms-, (24) 



that substitute now equations (1). Expressions (24) are 
valid at Ush ^ Uw, and under the conditions Rsh > Ro 



/M « 1{M,j/Mq) PC t> to - AO{Mej/MQ) 



3/2 



yr when the mass of swept-up gas is relatively large. 

Equations (14), (15) can be used for the estimates of 
maximum particle momentum when the instability is bal- 
anced by the wave advection and nonlinear wave dissipa- 
tion of Kolmogorov type. The asymptotic expression for 
Pmax(i) at high shock velocity subject to conditions (16) 



is the following (compare with equation (17)): 
Pmax _ 2M:aCcria)^creul^M'^/'^ ^ 



mc 



3 1/2 



8.9x10=^ (^,,/0.5)4/^f^ 



1/6 



tZ'- (25) 



mc 



The asymptotic expression at low shock velocity sub- 
ject to conditions (18) is 

Pmax ^ 0.2A>.a^ClMCUlrM''^<n ^ 4 g x 10^ 
^cr \ ^ ^7/3 / U«.,6 \ ^'^ ( B^r^ y 7/3 

0:5 j [jrj ivs^^^y "^^^ ' ^ ^ 

sX 9 n/2 (compare with equation (19)). 

In actuality the regime (26) can hardly be achieved 
since the size of the wind of a red giant star usually does 
not exceed 10 pc and the supernova shock breaks out of 
the wind cavity at t < 1 Kyr when the shock velocity is 
still high. Also, one has to keep in mind that the wind 
parameters are significantly changing over the course of 
star evolution. 

Equations (25), (26) show that, at the accepted values 
of parameters, the supernova shock moving in the wind 
of a red supergiant star can accelerate cosmic rays up to 
somewhat higher energies and the decrease of Pmax with 
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time is weaker than in the case of a shock moving in ho- 
mogeneous interstellar medium. Also, at a given age of 
supernova remnant, the value of Pmax is essentially higher 
for the explosion in wind. 



4. Discussion and conclusion 

The generation of turbulence by cosmic-ray streaming in- 
stability ahead of the shock is an integral part of cos- 
mic ray acceleration in supernova remnants. Wc consid- 
ered this process without limitation on the amplitude of 
random magnetic field and with allowance for the wave 
dissipation through Kolmogorov type nonlinearity and by 
ion-neutral collisions. The amplification of random mag- 
netic field at high shock velocity, Ush ^ 10'^ km s~^, leads 
to the increase of maximum energy of accelerated particles 
above the Bohm limit calculated at the interstellar mag- 
netic field strength. This is in agreement with the results 
of Bell & Lucek (2001 ). On the other hand, in support and 
development of earlier investigations of Volk et al. (J^.QSSJ , 
Fedorenko H1990|l . and Drury et al. H1996|l who dealt with 
not very strong turbulence and considered different non- 
linear processes, wc found that the nonlinear and linear 
wave damping may considerably suppress the level of tur- 
bulence at Ugh ^ lO'^ km s^^ and thus may decrease the 
maximum particle energy much below the Bohm limit. 
On the whole, the maximum energy of accelerated parti- 
cles with charge Ze rapidly decreases from > eV 
to much smaller energies ~ 10^"Z eV during the Sedov 
stage of supernova remnant evolution at the age from a 
few hundred years to about 10^ yr. Our main results are 
illustrated by Figures 1 and 2 for the cases of supernova 
bursts with kinetic energy of ejecta £ — 10^^ erg in hot 
rarefied interstellar gas (the gas temperature T = 7 x 10^ 
K, the average gas density including clouds uq = 10^^ 
cm"-^, the density of diffuse gas n = 3 x lO^'^ cm^'^, the 
interstellar magnetic field Bq = 2 fiG) and in the warm 
weakly ionized interstellar gas (T = 8 x 10'^ K, uq = 0.4 
cm^-^, n — 0.1 cm^-^, Bq — 5 fiG) respectively. The as- 
sumed efficiency of cosmic ray acceleration is ^cr = 0.5. 

It is remarkable that the highest energy -Emax which 
the cosmic ray particle can reach in the process of accel- 
eration does not explicitly depend on the value of inter- 
stellar magnetic field, see equations (17), (21). It is easy 
to check that the highest energy is attained at the time 
of transition from the free expansion to Scdov stage of su- 
pernova remnant evolution. Using equation (19), one can 
make the following estimate at t — to (Sedov solution is 
not applicable yet at this moment): 



2xaCcr{o)icrZmcu'^Ri) 

VaVgO 



8 X lO^^Z (uo/3 X lO-^km s^^)^ Mll^n^/^eY{27) 

at ^cr = 0.5, >f = 0.1, a = 0.3, no = n, the mass of ejecta 
Mej is in solar masses. The estimate (27) is in agreement 



with the estimates made by Bell & Lucek (2001). The 
corresponding random magnetic field comprises 

2aCcr{a)£, 

O-Dmax ~ — -t>0 ~ 

"J ya 

3 X 10~^^cr (uo/3 X lO^km s"^) n^^^G (28) 

in the upstream region of the size I?(pmax)/uo ~ ^Ro ~ 
O.li?o- The field can reach more than ~ 10^'^ G in the 
downstream compressed region behind the shock. The 
value of Smax given by equation (27) is by a factor of 
(5i3inax/-So ~ 100 larger than the Bohm value calculated 
for the interstellar field Bo- 

Equation (25) refers to the case of supernova explosion 
in the wind of progenitor star with the characteristic mass 
loss rate M = 10~^ Mq yr~^ and the wind velocity Uw = 
10®Mu,,6 cm s~^. 

The estimate of E'max in the case of supernova explo- 
sion in the wind of a progenitor star is: 



2xaCcria)^crZeulf^M'^^'^ 



1/2 



7 X lO^^Z (uo/3 X 10'*km s"^)^ y^M^s/u^^geV, (29) 

that in fact coincides with the estimate (27). In both cases 
the highest particle energy can hardly exceed lO^^Z eV 
even for the extreme value of shock velocity uq = 3 x 10*^ 
km s~^. 

Two essential assumptions are crucial for our conclu- 
sions. The first is the high efficiency of shock acceleration 
that transfers the fraction ^cr ~ 0.5 of the upstream mo- 
mentum flux entering the shock to the cosmic ray pressure 
and strongly modifies the shock profflc. This assumption 
is based on the results of numerical simulations of diffusive 
shock acceleration by Berezhko et al. I)199f)|l . Berezhko & 
Ellison 119995 ■ The second essential assumption is the high 
rate of nonlinear Kolmogorov-type dissipation of magne- 
tohydrodynamic turbulence that balance the cosmic ray 
streaming instability. This assumption is supported by the 
analytical and numerical investigations of MHD turbu- 
lence, see e.g. Goldreich & Sridhar l| 1995(1 . Verma et al. 
II1996II , and by the observations of Kolmogorov type spec- 
tra of turbulence in the interplanetary (Burlaga, I1995|l 
and interstellar (Armstrong et al.. I1995|l space. It should 
be pointed out that our description of nonlinear wave in- 
teractions is very simplified and can be considered as not 
more than an order of magnitude estimate. In particular, 
we ignore the complicated anisotropic structure of mag- 
netohydrodynamic turbulence in k-space. 

In the case of supernova remnant expansion into partly 
ionized interstellar gas, the wave dissipation by ion-neutral 
collisions was also taken into account. CoUisional damping 
may suppresses particle acceleration more than the non- 
linear wave dissipation does if the surrounding interstellar 
gas is not vastly ionized by the supernova burst (the case 
of Type la supernova) . 

The fast decrease of maximum energy of cosmic rays 
in supernova remnants with time may help to explain why 
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the shell remnants at the age more than a few thousand 
years are not prominent sources of TeV gamma rays pro- 
duced by very high energy cosmic rays (see e.g. Buckley 
et al. I1998|l in spite of the optimistic preliminary theo- 
retical estimates (Drury et al. I1994L Naito & Takahara 
I1994|l and in spite of a few detections of gamma-rays from 
these sources at about 100 MeV (Esposito et al. I1996|) . 
On the other hand, a number of young supernova rem- 
nants (SN1006, Cas-A, Tycho, Kepler, RXJ1713. 7-3946) 
were reported as cosmic ray accelerators of very high en- 
ergy particles with energies up to > 10 TeV based on 
the observations of nonthermal synchrotron X-rays or/and 
TeV gamma rays. (The effect for the old remnants is ev- 
idently compounded by the specific selection of the ob- 
jects near massive gas clouds that increases the genera- 
tion of secondary gamma rays through tt'^ decays whereas 
the high gas density no reduces Pmax as it is clear from 
equations (19), (22).) The variety of theoretical construc- 
tions by other authors were suggested to resolve this 
problem (Gaisser et al. [1998, Baring et al. [T^M Kirk 
et al. IWm Volk IMiTl Malkov et al. 2002") but the ulti- 
mate choice will be probably made when the data from a 
new generation of ground-based (CANGAROO-III, HESS, 
VERITAS, MAGIC) and space (GLAST) gamma-ray ex- 
periments in conjunction with new X-ray satellites will be 
available. 

The accounting for non-linear effects in the instability 
that accompanies the acceleration of cosmic rays in super- 
nova remnants may simultaneously eliminate two difficul- 
ties of the modern cosmic ray astrophysics. It raises the 
maximum energy of accelerated particles in young super- 
nova remnants above the standard Bohm limit through 
the production of strong random magnetic fields and thus 
helps to explain the origin of galactic cosmic rays with en- 
ergies up to ~ lO^^Z eV, though it does not cover energies 
larger than about 3 x 10^* eV even for Iron nuclei and this 
implies the extragalactic origin for cosmic rays with the 
highest observed energies. It also decreases the maximum 
energy of particles in late Sedov stage of supernova rem- 
nant evolution that allows to explain why these objects 
are not bright in very high energy gamma rays. 
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